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Abstract: Small fluctuations around a constant electric or constant magnetic field F 
are analyzed in a theory with pseudo scalar cf) with a coupling gtpFF. It is found that 
a magnetic external field leads to mass generation for the small perturbations, while an 
electric field suffers from a tachyonic mass generation in the case in which the field strength 
is higher than a critical value (related to the pseudo scalar mass). The vacuum energy can 
be exactly evaluated and it is found that an imaginary part is present when the external 
electric field exceeds its critical value. 
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1. Introduction 

1.1 Preliminary considerations 

A dynamical mechanism providing mass to vector gauge bosons is instrumental to match 
theoretical models of fundamental interactions with the particle spectrum observed in high 
energy experiments. A blueprint of dynamical mass generation is given by the Schwinger 
Model, or QED2, where fermions quantum fluctuations induce a mass term for the two- 
dimensional photon. Extension of this non-perturbative quantum effect to four-dimensional 
gauge theories has still to come because the gauge field effective action cannot be computed 
exactly in AD. In the meanwhile, the archetypal mechanism for gauge field theory mass 
generation is Spontaneous Symmetry Breaking, induced either by classical tachyonic mass 
terms [1] or by quantum radiative corrections [2]. The Coleman- Weinberg breaking of gauge 
symmetry avoids classical tachyonic mass terms and gives raise to a non-vanishing vacuum 
expectation value for massless scalar fields through radiative quantum corrections. In this 
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paper we are going to discuss a "complementary" mechanism, where mass follows from the 
breaking of rotational invariance induced by a classical background configuration of the 
gauge field strength. A real, or tachyonic, mass is obtained according with the magnetic, 
or electric, nature of the background field. The model implementing this effect consists of 
a scalar field <p non-minimally coupled to a U(l) gauge vector (non-abelian extensions of 
this model are planned for future investigations) through an interaction term of the form: 



The interaction term C\ has a long history dating back to the celebrated ABJ anomaly and 
neutral pion electromagnetic decay [3]. Moreover by keeping fixed the form and varying 
the strength of the coupling constant, C\ is equally well suited to describe the axion field 
currently appearing in many astrophysical and quantum field theoretical problems [4]. 

In what follows we are going to analyze the case in which the electromagnetic field is a 
purely electric or purely magnetic background, with special interest about the dynamics of 
its fluctuations. Before embarking this program and before giving a more detailed account 
of the main aspects of our approach, it is worth to recall some important steps already taken 
in the past in this direction. In the main part of this work we will stress more carefully 
analogies, as well as differences, with what we are proposing in this paper. In particular 
the fact that an external magnetic field modifies the dispersion relation of photons coupled 
to (pseudo)scalars was already discussed, for example, in [5]; there the authors have in 
mind an experimental set-up for the detection of pseudoscalars coupled to two photons 
based on the fact that the photon effective mass provided by the pseudoscalar coupling 
is responsible for an ellipticity in an initially linearly polarized beam 1 . Concerning the 
situation in which a background electric field is present, recently this problem has been 
analyzed in [8], where the authors show under which conditions an external electric field 
decays to pseudoscalars and discuss some particular configurations in which their results 
can be applied. Postponing a deeper analysis to what follows, we think that an important 
point to be stressed already at this early stage, is the fact that in the above studies the 
discussion is perturbative whereas, in the present paper, we are going to analyze a second 
order effective approximation for the dynamics of the fluctuations of the electromagnetic 
field, only after a full, non-perturbative treatment of the pseudoscalar. 

Before developing this part, we will shortly present some interesting features of the 
model in a naive form. We remember that, indeed, it is a special feature of the coupling 
(1.1) to generate physical masses, or tachyonic instabilities, when an appropriate classical 
configuration of the scalar field (j) is turned on: performing an integration by parts in the 
action associated with (1.1), we end up with an interaction term of the form 



We can then consider a background configuration selecting a preferred spacelike direction, 

e.g. (4>) = const. -S^x^, so that (1.2) gives a (2 + l)-dimensional, Chern-Simons type, mass 

1 This set-up is now a fully working project, PVLAS, at the I. N.F.N. Legnaro laboratories, taking data 
since May 1999: about this we refer the interested reader to [6] and references therein, as well as to [7] for 
a more direct connection with axion experiments. 



(1.1) 



(1.2) 



- 2 - 



term [9]. The resulting massive Chern-Simons model is embedded into a (3+l)-dimensional 
theory. Thus, rotational invariance in embedding space is broken [10]. It is interesting to 
compare this result with the case in which the background field is time dependent only, i.e. 
{(f>) = <p(t). The resulting Chern-Simons model is endowed with a tachyonic mass term 
[11], meaning that this type of background field is unstable. This kind of instabilities could 
play a role in baryogenesis, as it has been argued in [11], [12]. Moreover, a time dependent 
axion field may also affect the growth of primordial magnetic fields [13], produce effective 
Lorentz and parity violating modifications of electrodynamics and affect the polarization 
of radiation coming from distant galaxies [14]. 

1.2 A "naive" demonstrative approach 

In this paper we want to explore "the other side of the Moon", that is what happens if 
the gauge field strength F^ u , and not the scalar/axion field (p, acquires a non- vanishing 
background value. As an example, suppose (F^ u ) = B5^ 2 ^i/]3j i- e - (F^u) is a constant 
magnetic field along the x 3 direction. This choice is not ad hoc as it could appear: we 
are trying to mimic the QCD vacuum, where a constant, color, magnetic field lowers the 
energy density with respect to the perturbative Fock vacuum, where no gluons are present 
[15]. By introducing the proper kinetic term for (j> and the fluctuation field f^ u , we see that 
(1.2) turns into 

C = -\U»F U + ^ W + |50e 23 ^/ Q/3 . (1.3) 

Equation (1.3) shows as <fi couples only to the f 01 component of the fluctuation field 
strength. Accordingly, we can write the following effective Lagrangian in the (0 — l)-plane: 

An = -Ifoif 1 + ^ W + 9 -B4>f\ (1.4) 

By freezing the x 2 , x 3 dependence of the fields in (1.4) we get the bosonized Schwinger 
model [16] and solving the classical field equation for (f> we obtain the non-local effective 
action for the dimensionally reduced, massive, electromagnetic field 

An = -^/oi/ 01 - ^32 _ ^° 1 Q^ 01 ' (- 1 - 5 ) 

where the generated mass square is proportional to g 2 B 2 . Of course the full spectrum of 
excitations will contain more general types of spacetime dependence, so that we expect 
here an anisotropic mass generation effect. 

We can repeat the same calculations with an electric background, e.g. 

(Fpu) = ES\p,\oSi\ v ]. 

Then, instead of (1.5), we obtain for the small excitations in the (2 — 3)-plane 

C = "^23/ 23 + ^/ 23 0^3- (1-6) 

By comparing (1.6) with (1.5), we see that in the electric case there is a different sign in 
front of the non-local term, i.e. the generated mass is tachyonic. This situation for the 
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electric background can be partially cured providing a non vanishing mass to a to the axion 
field, so that the tachyonic mass generation sets in only around the critical electric field 
gE ~ m A . 

As we will see in the next sections, the theory can be solved beyond the two dimensional 
truncated sectors mentioned here. Furthermore, it is possible to compute the vacuum en- 
ergy resulting from small fluctuations around the chosen background. The vacuum energy 
develops an imaginary part for electric fields bigger than some critical value, of the order 
of mx/g. It can be worth now to remark a difference with respect to the instability in two 
dimensional QED. Both effects are triggered by an external electric field, but the standard 
Schwinger mechanism operates for any value of the constant electric field, provided it ex- 
tends to sufficiently large distances. Even with a weak field, one can perform enough work 
to create a new pair. By contrast, in our case, the effect is present only for field strengths 
bigger than a critical one, regardless of whether it extends over a large region of space or 
not. 

The paper is organized as follows. In Subsection 2.1 we set up the "sum over histories" 
formulation for the system under study: in particular we show how an effective Lagrangian 
for small fluctuations around a pure electric/magnetic background can be obtained, as a 
quadratic approximation, after integrating out exactly (and non-perturbatively) the scalar 
field and neglecting higher order terms. The next Subsection (2.2) deals in more detail 
with the mass generation effect, fully exploiting the difference between the electric and 
the magnetic case and determining the propagator and its spectrum in momentum space. 
Then in Section 3 the expression of the vacuum energy is introduced together with the 
exact expression for the propagator. The main steps in the computation of the free energy 
are outlined in Section 4 and the emergence of an imaginary part in the pure electric 
case is analyzed. Discussion and conclusions are in Section 5, with particular emphasis 
about the relation between the tachyonic mass-generation effect and a possible instability of 
homogeneous electric fields above some threshold. Three appendices follow, where technical 
details can be found about the solution of the eigenvalue equation for the momentum space 
propagator (Appendix A), the determination of the propagator itself (Appendix B) and 
the explicit computation of the free energy (Appendix C), with particular attention at the 
behavior in the infrared and ultraviolet limits. 

2. General analysis of the small perturbations 
2.1 Path integral 

An effective technical framework to investigate quantum fluctuations around an external 
background configuration is provided by the Feynman path integral formalism. The parti- 
tion function(al) encoding the dynamics of the interacting (f> and fields reads 




(2.1) 



where 



C = -\f, v + 9 ^e^F, v F pa + \d^d»<j> + 



(2.2) 
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and both the gauge fixing and ghost terms are, momentarily, understood in the functional 
measure [DA] . For the sake of generality, we assigned a non- vanishing mass to the pseudo 
scalar field <j>. Since the path integral (2.1) is gaussian in 0, the scalar field can be integrated 
away exactly, i.e. 



J [£>0]exp j-2 J d 4 x 

'U + m\ 



det 



-V2 r „2 



expj^ | d 4 xe^F^F pCT ^-^ e ^F Q/3 F 7<5 |, (2.3) 



where ji is a mass scale coming from the definition of the measure \P<j>\- Integrating out 
the (f) field induces a non-local effective action for the A field. We stress that this effective 
action, being obtained at the non-perturbative level, takes into account (as an effective 
action for A) all the effects due to the presence of the pseudoscalars at all perturbative 
orders. This is a crucial difference with many of the previous works on the subject. As 
is clear from (2.3), the resulting path integral is quartic in A but, even if it cannot be 
computed in a closed form, the background field method provides a reliable approximation 
scheme to deal with this problem. We thus split F^ u in the sum of a classical background 
(F pi/ ) and a small fluctuation f pu : 

F^ = (F^) + V (2.4) 

In the case of a pure electric or a pure magnetic background we have 

e^ a(S (F flu )(F a( ,) = 0. (2.5) 

Thus, expanding C up to quadratic terms and dropping a total divergence, we obtain for 
the effective Lagrangian of the fluctuations of the electromagnetic field 

C {2) = -\(F, u )(Fn - \f»„r - f^^(^) e ^(F pCT )/^^_/ 7 ,. (2.6) 

As we will see in a while, there is an important sign difference between the magnetic and 
the electric case. 

2.2 The signature of the mass generation: ordinary versus tachyonic 
In this subsection we turn to the analysis of the contribution 

e^(F^y^ s (F pa ) (2.7) 

of equation (2.6), which will be responsible for ordinary or tachyonic mass generation 
according to whether (F pu ) represents an external magnetic or electric field. We notice 
that our results are not inconsistent with those that can be deduced from equation (31) of 
[8] , from which it is clear that the sign of the contribution from the interaction Lagrangian 
changes in the case of purely electric or purely magnetic background. 
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Let us now start from the case of a constant magnetic field B. Then, without loosing 
generality, we can rotate the reference frame to align an axis, say x 1 , with B. Accordingly, 
(F^) = B5 M2 K]3, so that 

e puaf3 (F pi/ )e p ^ s (F pa ) =4B 2 e 23a/3 e 23 ^ 5 . (2.8) 

It is clear thus that a, (3, 7 and 5 in (2.8) can take only the values 0, 1. Denoting by O 2 )^"^ 
the 2x2 Minkowski tensor, we must then have that 

£ a/323 e7 523 = A ^(2)^(2)^8 _ (2)^8(2)^ 

= ^fe<o)-<o)^ 7 o))> (2-9) 

where P^ is the projector onto the (0 — l)-plane, i.e. P^ = Wrf* 1 . Contracting in 
relation (2.9) the couples of indices (° 7 ) and ( l3S ) we obtain 

e a/323 e af 3 23 =A(2x2-2) = 2A; 

since in lowering the indices a/3 there is the time involved, e a ^ 23 e a /3 23 = ~~ 2, so that A = — 1. 
Analogously in the electric case we can take F^ u = E5[^ 5i^. Then 

e a(3fMP '{F pv )e^ pa (F pa ) = 4E 2 e afi01 e rm , (2.10) 

with the indices a, (3, 7, 5 taking only the spatial values 2, 3. Now, following the same 
procedure and observing that no minus signs are involved in lowering spatial indices, we 
find 

e a/301 e 7<501 = ^ (2) gory (2)^8 _ (2)^a5(2)^ 7 ^ 
~ ^(23)^(23) ^(23)^(23)' 

again we introduced the projector notation = ^5^ v , where ^ 5^ u is the 2x2 Kronecker 
delta. 

We can now write the equation of motion for the fluctuations, which from the action 
(2.6) turns out to be 

V" = 8g 2 e^(F pu )— * e P^ {F R/ (2 . n) 

□ + m\ 

The above equation can then be expressed, for both the electric and magnetic case, in 
Fourier space and in terms of the vector potential of the fluctuations, which we will call 
a p . Since f^ v = d p a v - d u a pi we have 



ka (k^a u - k u a^) = K 9 (k 2 g pu - k"k u ) o„. (2.12) 
k 2 — m A 



Here and in what follows we define 



& = Pfik a = g* a k a (2.14) 
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and 

in the magnetic case, or 

in the electric case. Equation (2.12) can be written also as 

where we have defined 

V~ llxv {k) = (k 2 g^ - k^k u ) - — ^ — ~~ — 2~ (Pg^ - k^k v ) . (2.15) 

The operator in (2.15) can be diagonalized 2 as shown in appendix A. Four linearly inde- 
pendent physical states satisfy the eigenvector equation of V~ ltiu (k); we will call them k^, 
k^, k~h, E^. Their definitions and corresponding eigenvalues are enlisted in table 1. There 



Eigenvalues 


Eigenvectors 





k^ 


j 2 hik 
k 2 — m\ 


( e^ am k a in the electric case 

h = \ 

[, e 23 ^ a k a in the magnetic case 


k 2 




k 2 


E n = e wk u k p k± 



Table 1: Eigenvalues and eigenvectors of V ltiu {k). 



are, thus, two non-zero eigenvalues and one of them, k 2 , has degeneracy two; moreover the 
"gauge" eigenvector is associated to the zero eigenvalue. This last results changes if we 
study the problem in the covariant a- gauge, when the eigenvector k^ is then associated to 
a non-zero (but still background independent) eigenvalue, 1/q (the corresponding inverse 
propagator will be called V^ 1 ^ u {k;a) as defined in what follows). 

As promised above, we can now compare the obtained results with those derived in 
previous works on the subject. Again, we observe that a key point in our derivation is 
the (non-perturbative) path-integral procedure used to obtain equation (2.3): in this way 
all the effects due to the quantum fluctuations of the pseudoscalar field, at all orders, are 
taken into account. In particular, if we concentrate on the purely magnetic case, it is 
then clear that this constitutes a generalization of the results obtained in [5], where the 
secular equation is obtained considering plane wave solutions to the classical equations of 

2 We will summarize here the final result about eigenvectors and the corresponding eigenvalues with all 
the relevant definitions, referring the reader to the mentioned appendix for the detailed computation. 
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motion associated to the Maxwell + Klein-Gordon action for the coupled electromagnetic 
and pseudoscalar fields. An analogous result for the electric case in presence of massless 
pseudoscalars can, for instance, be found in [17]: here the dispersion relations are obtained 
under physically very sensible restrictions but are, anyway, of perturbative character. Also 
the more detailed analysis of [18] uses a different kind of approximation with respect to 
the one employed in our calculation since the starting Lagrangian in equation (1) of [18] 
is different from (2.6), our quadratic approximation to the full, non-perturbative, effective 
result in equation (2.3). We can thus trace back the differences between our results for the 
eigenvalues of the propagator and the one already derived in the literature on the subject, 
to the fact that we have taken into account the effects due to the pseduscalar fields in a 
substantially non-perturbative way. 

3. Path integral quantization of small perturbations: propagator and vac- 
uum energy 

To compute the vacuum energy we proceed further in the covariantly quantized fashion we 
started in the previous section. In a covariant a-gauge the path integral for the partition 
function can be rewritten, using (2.15) and going to Euclidean space, as 

Z = [det (□ + mi)] - 1/2 J [Vf] J [VA^S [d^ - /] e"£ / ^ x det [d 2 ] x 

'4 1 



m A 



A v (k) } . (3.1) 



Here we have already performed the functional integration over 0, we remember that we 
are approximating the higher order Lagrangian up to the second order in the fluctuations 
and understand the path-integrals in the Euclidean sector. The functional integration over 
/ can now be done and we get 

det [d 2 ] J[VA^ V {- J ^A,{k)V-^ v (k- ia )A v {k)\ 

Z = L — m L , (3.2) 

[det (D + mi)] 1/2 

where 

V- llM/ (k;a) = (k 2 g^ - Wk u ) - , _ " 9 (k 2 g^ - k"k u ) + -WW = V^ v {k) + -Wk u . 

' k 2 — m z A 'a a 

In order to solve (3.2) we have to find the eigenvalues, of V~ 1 ^ iu (k; a): Z equals then the 
product of these eigenvalues. The requested eigenvalues are those related to the physical 
polarizations found in section 2.2, i.e. k^, kj[, E^, with eigenvalues k 2 — nk 2 /(k 2 — mi), k 2 
and k 2 respectively. As already discussed at the end of the previous section, the eigenvector 
W has now a non-zero eigenvalue 1/a. 

The vacuum energy, or free energy, is then 

3 
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so that 



+Y fJw> ihl{k2 ~ mlh (3 - 3) 



The propagator, i.e. the inverse of X> 1/liy (/c;a), can also be found exactly, as shown in 
appendix B, and results to be 

[k- a) - k2 [9 k2 ) + k2 ^ 2 ^ fc2 _ m 2 ) _ ^ + 

"77^- (3-4) 



fc 2 (fc 2 (fc 2 -m^) - K fe 2 ) (fc 2 ) 2 ' 

4. Evaluation of the vacuum energy and of its imaginary part 

To compute the vacuum energy, apart from more standard contributions, we have to com- 
pute the following integral 

elk -, /, 2 fck 



UK k; k, m A ) = J j^y[ In \k - j . (4.1) 

It can be evaluated in closed form and we first rewrite it as 

/d 4 k 
-^—^ In (k 2 - m 2 A ) , (4.2) 

so that we can separate the second common contribution, from the first one, i.e. 

I(k, k; k, m A ) = J -0^ In [k 2 (k 2 - m 2 A ) - K~k 2 ] . (4.3) 

Using (4.2) for the right hand side of (4.1) appearing in (3.3), we get for the free energy 

VT (' d^h 

W = V —I{K k k, m A ) + VT J m (k 2 ) ■ (4.4) 

The computation of I(k,k; k, m\) is performed in appendix C. The infrared behavior 
is then extracted in appendix C.l. Ultraviolet divergences are regularized by a cut-off A and 
the leading contributions to the free energy are computed in appendix C.2. Thus the final 
result for the vacuum energy density is obtained multiplying (C.25) by l/(47r 2 ) (because of 
(C.l)) and substituting for I(k,k; k, to a ) in (4.4), where the second contribution can also 
be exactly evaluated. The final result is 



VT 
W = 

8vr 2 L 



/i 4) A 4 + tfU 4 In A + /i 2) A 2 + lP\ A 2 + J<°> + iffi 



(4.5) 
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where the various contributions are defined at the end of appendix C. We will be especially 
interested in VTI^ /(8ir 2 ), which, according to equation (C.26), is 



/( o) = YL 

8vr 2 



49k 2 + 132rem 2 A + 132ml re 2 + 3rem 2 A + 3mi , 
a a H a — In 2+ 

576 48 
-^7^1n(K + mi) + ^-lnK) 



(4.6) 



From this expression we see that the vacuum energy acquires an imaginary part in the case 
K + m\ < 0, i.e. when the tachyonic modes are present. If we make use of the prescription 
re + m\ — ► re -\- m\ + ze, the value of the imaginary part is 

Notice that, as opposed to the real part, the imaginary part of the vacuum energy is 
cut-off independent. This is so because only the infrared region of the integrand in (4.3) 
contributes to the imaginary part. 

We also observe that all contributions in (4.5) are finite in the limit of vanishing external 
field, i.e. when re — ► 0. This is immediately evident for the terms lf \ An A an< ^ 

IfaAi from their expressions at the end of appendix C. Concerning the contribution 
the two divergent terms of opposite sign in the second line of (4.6) give a finite contribution 
in the limit, as shown in (C.12). 

The cut-off dependent parts are regularization dependent, as it is known from general 
experience with the regularization of divergent integrals. In particular subleading diver- 
gences (i.e. the A 2 terms in (4.5)) are highly dependent upon the regularization scheme and 
may vanish in certain of them. In the present problem, due to the complexity of the inte- 
grand in (4.3), other regularization schemes, like £- function regularization or dimensional 
regularization, are difficult to implement. 



5. Discussion and conclusions 



In this paper we have discussed effects of mass generation in an external magnetic field 
and of tachyonic mass generation in an external electric field in the case in which there 
is a pseudo scalar field with a pseudo scalar coupling g(f)F fll/ F flu . The effects due to the 
presence of the psedo scalars are fully considered: indeed in the derivation of an effective 
theory for the fluctuations of the electromagnetic field, the contributions from quantum 
fluctuations of the pseudo scalars are taken into account using a non-perturbative approach. 

In the purely magnetic case the mass production for the electromagnetic fluctuations 
can be interpreted as in the cited works already present in the literature: the differences 
that we find in the eigenvalues of the propagator can be traced back to the non-perturbative 
character of our approach, as opposed to the perturbative analysis performed elsewhere. 

A more careful discussion is instead required for the purely electric background, in 
connection with what we have called tachyonic mass generation. Indeed the appearance 
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of an imaginary part in the free energy suggests the presence of an instability for homo- 
geneous electric fields beyond some threshold due to pseudo scalar coupling. This means 
that the vacuum state must be redefined, to obtain a correct ground state for the theory. 
The analysis of what becomes the true ground state, a genuinely non perturbative effect, 
is beyond the scope of this manuscript. Nevertheless we would like, at least, to suggest a 
possible and simple, although incomplete, answer to this question. Indeed an analysis of 
Cornwall [22] indicates that in the case of a 3-dimensional, Euclidean, tachyonic mass term 
(which in our discussion is generated by the time dependent scalar field of equation (1.2) 
but can have also another origin) the final result is the formation of an inhomogeneous 
state. Then, properly generalizing to our set-up the results described in [8], where inhomo- 
geneous electromagnetic fields are shown to decay in pseudoscalars, it is not unreasonable 
to understand under what we have called tachyonic mass generation a real "tachyonic 
instability of the vacuum". 

From the result of equation (4.6) it is possible to see that the free energy is well-defined 
in the limit of vanishing background fields (k — ► 0). Thus our effect is a genuinely non- 
perturbative one and does not relate to a particular choice of the regularization scheme. 
Moreover it is worth pointing out again that this "tachyonic instability of the vacuum" 
for fluctuations around a constant external electric field, is characterized by a threshold 
effect, i.e. the tachyonic mass generation is switched on for electric fields high enough, so 
that k + m\ < 0. In the case of the neutral pion, we can obtain the value of the effective 
pion-photon coupling, defined by equation (1.2), from the observed value of the neutral 
pion lifetime [20] and from the value of the decay rate given the coupling (1.2) [21]. This 
gives us the values 

g = 2.53 • 1(T 5 MeV -1 , 
m n = 134.97 MeV, 
E crit . = (1 GeV) 2 . (5.1) 

This is a very high electric field, not available in normal laboratory conditions. Further- 
more, if it were available, it would reveal the composite structure of the pion and the 
effective g^F^F^ coupling, used here, would not be applicable any more. 

A different question would be then the study of this effect in the case of hypothetical 
axion particles. In this case, the threshold for the tachyonic mass generation to be set up 
becomes lower as lower values for the mass of the axion are considered. 

Apart from the purely electric case, which is more subtle and, maybe, more exciting 
because of the exotic tachyonic mass term, it could be interesting a more detailed analysis of 
the magnetic case in connection with the set-up of PVLAS a presently running experiment 
at the Legnaro I. N.F.N, laboratories, near Venice, Italy. 

Finally a totally different role for these effects, could be in the context of QCD. There, 
it is known that an external chromo-magnetic fields presents tachyonic instability. If we 
were to add a particle with coupling to ^ ua ^ F^ V F^ (this particle could represent a pseudo 
scalar bound state of quark and anti-quark pairs), we know that the effect of the external 
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chromo-magnetic field together with the pseudo scalar coupling is of generating mass. The 
interplay of these two effects could then be an interesting subject for further research. 
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A. Eigenvectors and eigenvalues of V 1>iu (k) 

We will now find the solutions to the eigenvector equation 

V- liiV {k)a v = X k a^ (A.l) 

where V~ 1 ^ u (k) is defined in equation (2.15). 

Firstly, there is a "trivial" gauge solution = H(k)k^, since k^ is orthogonal to both 
terms enclosed in round brackets in the definition (2.15) of T>~ 1 ; this can be seen from the 
relations 

(k 2 g^ - k»k v ) fc„ = 

and 

(k 2 g» u - k»k v ) k^ = 

using the properties g^ u k v = g^ v k v and k v k v = k v k v . 

A second, non trivial polarization is k^ = e^ a01 k a in the electric case {k^ = e 23 ^ a k a in 
the magnetic case), since k^k^ = = k^k^ thanks to equations (2.13) and (2.14); a M = k^ 
is associated to the eigenvalue 

k 2 - -J^. (A.2) 

fc z — m\ 

Then, remaining independent physically relevant polarizations must be orthogonal to 
both k^ and k^, hence they can be parametrized as 

a M = ^"^djk a kp; (A.3) 

since given by (A.3) is not affected by the "gauge transformation" 

d v — ► d u + Xik u + X 2 k u , 

it then follows that only two components of d v are physically relevant. One of these, 
which we will call k-^, is obtained when d v = k v : it is thus orthogonal to k^ 1 , k^, k^ and 
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corresponds to the eigenvalue k 2 . This is the same eigenvalue of the last eigenvector, E^ 1 , 
which is given by 

This can be verified by inserting it into (2.12) and observing that 

k 9 nv k^k v oc k^k v . 

The above equality holds because 

k 9 [iv ~ k^k v 

is a projector onto the space orthogonal to k^ in the 2- dimensional subspace, where the 
direction orthogonal to k^ is nothing but k u . 

B. Determination of the propagator 

To determine the propagator in the covariant a-gauge, T>~ liiu (k;a), with the four dimen- 
sional quantities that are at our disposal we consider the ansatz: 

A[g^- ^-f- + Bg^ + Cfrk" + D- 



k 2 J * {k 2 f 
To determine the coefficients A, B, C, D we can now compute T>~ lfMU (k; a)T>~ 1 l/a (k; a): 



1 



(k 2 g^ - k»k u ) - -5 T (k 2 t v ~ k^k u ) + -fcT 

k — mi a 



A [g 



ll/a W 



+ Bg va + Ck v k a + D 



a 
kijka. 



— Ak 6^ Ak^k a —r ~- (k k^k a j -\- k^k^ ~\~ 
k A — mi a 



A.k^ k(% ~t~ t4./l^/cq, ~\~ 



kA (k 



k 2 



m, 



A. 



+Bk 2 5£ - Bk^ka 



(k k^kf^) -\- k^k a -\- 

OL 



k 2 — m\ 



~\~C k k^ k(% C k k^ k^ 



+D- 



k^k r 



k 2 



D- 



k»k r . 



k 2 - 
kD 



C - 

(h 2 h^h —h 2 h^h \ 4- —h 2 h^h 4- 

I rv rv /Vq, rv rv /Vq, I iv rv Ajq, ~ 



a 



k 2 k 2 



m> 



, iv rvr\ 



a k 2 



Now comparing terms with the same tensorial character 



8£ : Ak = 1 

k^k a : — t- = A 
a k z 



A: 2 - 



k^k a : -B - Ck 2 = 



k 2 



m, 
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An independent subset of (four of) these gives (consistently with the remaining equation) 
the final result for the coefficients: 



A 
B 



C 
D 



1 

fc2 



Kk 2 A 



k 2 {k 2 — m\) — nk 2 



k 2 (k 2 (k 2 - m 2 A ) - nk 2 ) 

K 

k 2 (k 2 (k 2 — m\) — nk 2 ) 

Q, 



so that we can finally write 



Wk v 



+ 



Kk 2 g^ 



k 2 (k 2 (k 2 - m\) - Kk 2 ) 



+ 



Kk^k 1 



+ a- 



k^k h 



k 2 (k 2 (k 2 - m 2 A ) - Kk 2 ) (k 2 ) 2 ' 
C. Explicit computation of the vacuum energy integral 

We concentrate in this appendix on the computation of the contribution I(k, k; n, m) of 
equation (4.3). Of course since the integral is divergent, it must be properly regularized 
and we choose to do that by putting an infrared (e) and an ultraviolet (A) cutoff on the 
modulus of the momentum k and of its projection k, exploiting some of the arbitrariness 
in the choice of the regularization scheme. Thus integrals written with implicit integration 
domain, like f d 4 k(. . .), are to be understood as performed in the domain of the variables 
(ko, ki, kii ks), which is the inverse image of the domain e<r<A, 0<??< 2tt, e < p < A, 
< w < 2tt in the variables (r, p, w) under the following change of variables in Euclidean 
space: 



k 1 
k 2 
\k 3 / 



T(r,$,p, w) 



/ r cos \ 
r smfi 
pcosw 

\ p sin w ) 



with Jacobean 



JT 



/ cos-d 
sini? 


V 



— r sini? 

T COS 1? 









cosro 
sintu 






— psinzu 
pcosro . 



whose determinant is 



so that 



det (JT) = rp 



d 4 k = rpdrdpd'ddw . 
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We thus have 



I(k, k; re, toa) 



/ 



^ln[A: 2 (fc 2 -mi)-^ 2 ] 

2tt f -2TT f-A f-A 



Ztv rl\ pi\ 

<h) / dw Jdr J dp det (JT) x In [(r 2 + p 2 )(r 2 + p 2 - m\) - rep 2 ] 



(2vr) 4 Jo jo j £ 
1 f A f A 

= -p^ Jdr Jdprp x In [r 4 + (2p 2 - m 2 A )r 2 + p 2 (p 2 -m\- re)] 

We have thus the problem of computing J(r, p; K,m\), which can be calculated in closed 
form: the final result can be expressed as 

p 2 [re + 3(m 2 - 6r 2 - p 2 )] 
J(r,p;re,m A ) = — h 

-B 3//2 (C-s/B\ A 3 / 2 / k + C-Va\ 

\k 2 + 3re(mi - 2r 2 ) + 3(m A + 2D)1 , , 2 , 
+- — ^ — In (-rep 2 + D) 

if we set 

A = A(r; jtia, re) = (re + "J-a) 2 ~~ 4rer 2 

£? = B(p; niA, re) = 4rep 2 + m\ 

C = C(r, p; m A ) = mi - 2(r 2 + p 2 ) 

-D = D(r, p; m A ) = (r 2 + p 2 ) 2 - m\(r 2 + p 2 ) = + up 2 . 

The final stage consists in evaluating it in the infrared (e — > 0) and ultraviolet limits. In 
this case this is equivalent to the computation of 

lim J(r, p; re, m\) 

(r,p)-(0,0) 

and the extraction of the divergent contribution in 

J(r = A, p = A; re, to a ) 
when A — > oo. This is done in the next subsections. 
C.l Infrared limit 

In this subsection we consider the infrared limit and in the following we will use the symbol 
"=" to imply that two expression are equivalent in the infrared limit, i.e. they have the 
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same limit. Moreover we will get rid of the square roots by means of the following results 

2rer 2 



(re + m\) 2 
V mi 



1 - 



re + m 2 A \ 2 



1 + 



2re/? 2 



m. 



(C.3) 
(C4) 



We now turn to the contributions in the various lines of equation (C.2). The first one gives 
no problem: 

p 2 [re + 3(m 2 A - 6r 2 - p 2 )] 



24 



^ 0. 



(C5) 



The second one has a well behavior in the B term; more care has to be paid in the 
logarithm: 



B 3/2 

48re 



In 



C 



rrij: 
48k 



In 



m\ - 2(r 2 + p 2 ) - \/4kp 2 + m\ 



ra\ - 2(r 2 + p 2 ) + ^4np 2 + 



m, 



^ A In 



48re 



m\ - 2(r 2 + p 2 ) 



1 + 



2np z 



l-2{r 2 +p 2 )+ml (l + ^) 



m 



A 



18k In [-m 2 A (r 2 + p 2 ) - rep 2 ] - ^ In [(^i) 2 ] . (C.6) 

The last term on the second line of (C.2) again gives troubles only inside the logarithmic 
term, which can be elaborated as follows 



In 



re + C- y^' 
k + C + \/A, 



In 



In 



re + m\ - 2(r 2 + p 2 ) - y^re + m\) 2 - 4rer 2 
re + - 2(r 2 + p 2 ) + \J(k + m\) 2 - Anr 2 
re + — 2(r 2 + p 2 ) — |re + m^J ( 1 



__2Krf_ 

|k+to^J 2 



= < 



In 



In 



re + m 2 - 2(r 2 + p 2 ) + |re + m|| (l - j^r^) _ 

if re + m 2 ^ > 



_( r 2 + 2 )+ ^H^_ 



( K +m 2 )-(r 2 +p 2 )- 
/ ~(r 2 +p 2 )+- 







( K +m 2 )-(r 2 +p 2 )- 



ln 



-(r 2 +p 2 ) + 



. .2 



if re + m a < 



SignfK+m^)' 



(re + m|) - (r 2 + p 2 ) 



Since we also have 



Sign(re + m 2 A ){ln[-m 2 A (r 2 + p 2 ) - up 2 ] -ln[(re + m|) 2 ] } .(C.7) 
A 3/2 ^ (re + m^l 3 



48re 



48re 



(C.8) 
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the two previous results, (C.8) and (C.7), combine in a neat way: the sign in the first factor 
exactly combines with the absolute value of the second factor 

\k + m^| 3 Sign(ft + m A ), 

so that 

48k y K + c + ^fA) 

~ ln t" mi(r2 + p2) ~ ^ ~ ln t (K + mi)2 ] • (C - 9) 

For the last term in (C.2) we do not have too much work. The factor before the logarithm 
has no problems and we can simply forget about the r and p dependent parts. Instead 
inside the logarithm we can neglect higher order terms in the limit we are interested in, so 
that 

[k 2 + 3 K {m 2 A - 2r 2 ) + 3(m\ + 2D)] ^ 



ln (-k P 2 + D)^ 



48 

= ^ + ^§ + ^ I-™ \{r 2 + p 2 ) - np 2 ] . (CIO) 
The desired result, J(0,0;K,m A ) is then (C.5) + (C.6) - (C.9) + (CIO), i.e. 

J(0, 0; k, m A ) = ^ {In [-m 2 A (r 2 + p 2 ) - k P 2 } - ln [(mi) 2 ] } + 

- {K + ™} ? ln[-m 2 A (r 2 + p 2 )- K p 2 ] + 
+ (^+m|)3 ln[(K + mi)2] + 

K 2 + 3Kf77.il -\- 3l7T.l r 2/2 2\ 11 

+ ^ 4 In [-m 2 A (r 2 + p 2 ) - K p 2 } 

= (^ ln(B+mi) _Ml! lnW) . (cm 

As also pointed out in the main text this contribution is finite in the case of vanishing 
external field (k — > 0), since 



lim 

K->0 



( K + <> 3 1 „( K + ,„i)-M£lnK) 



24k v a/ 24k 

C.2 Ultraviolet limit 



4 

m 



[ ln(m 6 } + x] (C12) 



To tackle the problem of the ultraviolet behavior of the energy density, we analyze the 
limit in which r — > oo, /? — > oo. As already discussed we will first set r = p = A and 
then approximate the various quantities as A — > oo. For the relevant expressions, already 
encountered above, we get 3 

* = -*A'(l-^). A> -<!i±<£). ... «,») 



3 We now use the ~ symbol to imply the same behavior in the ultraviolet limit. 
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B = 4kA 2 1 + 



m A 
4kA 2 



B z ~ 16k 2 A 4 ( 1 + 



m> 



C = -4A 2 1 - 



4A 2 / ' 



2kA 2 / ' 



k + C = -4A 2 1 - 



4A 2 



(k + C) 



-l 



4A 2 



D = 4A 4 1 



2A 2 / ' 



log(-« + £>) ~ log(4A 4 ) - 
Moreover we also have the well known expansions 



k + m A 
4A 2 



4 k + 2m\ (k + 2m 2 , 



(C.14) 



(C.15) 



4A 2 



32A 2 



X X 

Arctanh (x) = x + — + — + 0(x 7 ) 

o 

arctan (x) = x — — + — + C(x 7 ), 



(C.16) 

(C.17) 
(C.18) 



which we are going to use in the following. In particular we can consider in generality the 
expansion of the following expression 



w 3/2 ( z — W X I 2 

m 



1/2 



+ W 1 / 2 



w 



= < 



= < 



3 / 2 Arctanh 

V z ) 

' M 3/2 Arctanh (H^l) if w > 

=W/ 2 = \ w \l/2 and w 3/2 = | w |3/2 

M 3 / 2 Arctanh if w < 

tt; 1 / 2 = l \ w \ 1 / 2 and ttJ 3 / 2 = — «|tt;| 3 / 2 
' M 3 / 2 Arctanh (H^l) if w > 

z|w| 3 / 2 i arctan ^^—7—^ if w < 
' |H 3/2 Arctanh (^) if > 

M 3/2 arctan (^1) if w < 
,3/2 ,H 1/2 



. , ,1 /k| 1/2 \ 
sign(w)- I — — I + 



+ 



1 / \w\ 1 ' 2 



sign(iw) 



111 



1 1/2' 



+ 



so that 



3/2 



In 



2 — w 



1/2 



z + u> 1 / 2 



0,00 



(2n + l)^ 2 ^ 1 



(C.19) 
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From the above result, if we identify 



w < — ► B, 
z < — > C 

and we consider an overall l/(24re) factor, using properly (C.14) and the first equation of 
(C.15), we get 

S 3/2 (c-Vb\ reA 2 re 2 + 3rem 2 + 6m 4 



48k 



In 



C + ^B 



+ 



72 



(C.20) 



for the second term in (C.2). In the same way, starting again from result (C.19), together 
with the identifications 

w < — ► A, 

Z < ► K + C 

and taking into account an overall — l/(24/c), (C.13) and (C.15), we obtain for the third 
term in (C.2) 



A 3 / 2 , f K +C-y/A 

■ In 



48k 



k + C + y r A 



kA 2 4k 2 + 9rem 2 + 6m 4 



6 



72 



(C.21) 



The last term in (C.2) has also a logarithmic part and, using (C.16), can be approximated 



as 



[re 2 + 3Km 2 A - 6rer 2 + 3m\ + 6D] In (-up 2 + D) ~ 

~ " 2 + 3 « m A + 3 ™A ln ( 4A 4) _ K + 2m A A 2 ln ( 4A 4^ + ln ( 4A 4) + 



48 

, ( K + 2mj) 2 _ re+jmj ^ 



64 



(C.22) 



We are now concerned with the easiest term in (C.2), namely the first, for which we get 



p 2 (re + 3 (m 2 A - Qr 2 - p 2 )) re + 3m\ 2 7 



24 24 
Summing up equations (C.20), (C.21), (C.22), (C.23), we obtain 



A 4 



A 4 



21nA + ln2 



7 



+ A 2 



2k 



+ 3m\ re + 2m 2 , 

^ ^ (2 ln A + ln 2) 

24 4 

2 \2 .„2 



+ 



(C.23) 



(C.24) 



5k 2 + 12nm\ + 12mA U + 2m\) 2 k 2 + 3/emi + 3mA , , , , „ N 

H + — H (2 ln A + ln 2) . 

72 64 48 

C.3 Result for I(r = A, p = A; re, ttja) — lim^p^^o) A r > P! K i m A) 

From the results (C.ll), (C.24) of the two previous subsections we obtain what we are 
interested in, 

J(r = A, p = A; re, m A ) — lim J(r, p; re, ttja) ~ 

0*,pH(o,o) 
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A 



2 In A + In 2 



+ A 5 



2k + 3m\ k + 2m 2 A 
24 4 

,2 \2 .„2 



(21nA + ln2) 



+ 



5k 2 + YlKm\ + 12mA (n + 2m\) 2 k 2 + 3Km\ + 3mi , , , . 

^ ^ + — H (2 In A + In 2) + 

72 64 48 v ; 



72 

(k + m^) 3 
24k" 



In (k 



+ m A ) + ^ln(m A ). 



~ (/i 4) " In 2 + |) A 4 + (4 4 A - 2)A 4 In A + if A 2 + /£ A A 2 + / (0) + 4° A (C.25) 
where for convenience in the last line we have used the following definitions: 



j(4) _ 
7 lnA — 4 



21n2-- 



r (2) 

1 A — 



r(2) _ _ 
i lnA — 

/(0) 



2k + 3m 2 , k + 2ml , 

7T, — : — - ln 2 

24 4 

k + 2m\ 



49k 2 + 132Kmi + 132m 4 , k 2 + 3nm\ + 3m\ , 
= — = ^ H -e In 2 + 



576 48 
(k + tti 2 ,) 3 , 2 , (m 2 A ) 3 



r(0) _ 
J lnA — 



24k 

k 2 + 3nm\ + 3m\ 
24 



24k 



ln (m 2 A ) 



(C.26) 
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